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Recently, several authors have investigated topological phenomena in periodically-driven systems 
of non-interacting particles. These phenomena are identified through analogies between the Floquet 
spectra of driven systems and the band structures of static Hamiltonians. Intriguingly, these works 
have revealed that the topological characterization of driven systems is richer than that of static 
systems. In particular, in driven systems in two dimensions (2D), robust chiral edge states can 
appear even though the Chern numbers of all the bulk Floquet bands are zero. Here we elucidate 
the crucial distinctions between static and driven 2D systems, and construct a new topological 
invariant that yields the correct edge state structure in the driven case. We provide formulations 
in both the time and frequency domains, which afford additional insight into the origins of the 
"anomalous" spectra which arise in driven systems. Possible realizations of these phenomena in 
solid state and cold atomic systems are discussed. 



The discovery of the quantized Hall effectP revealed 
the existence of a powerful new class of extremely robust 
quantum phenomena which can be observed with high 
fidelity, largely independent of sample size, shape, and 
composition, up to macroscopic (millimeter) scales. The 
robust nature of these phenomena can be linked to the 
presence of energy gaps for bulk excitations, combined 
with the existence of nontrivial topological structures as- 
sociated with the systems' ground state wave functions^. 
Recently, these ideas were used to predict the existence 
of new classes of materials^', the topological insulators, 
which were found experimentally^^ shortly thereafter. 

Meanwhile, a wide variety of new experimental tools 
has been developed for actively controlling and probing 
the behavior of electronic and cold-atomic systems. For 
cold atoms, various methods have been proposed for cre- 
ating "synthetic gauge fields," which mimic the effects of 
magnetic fields^*] r spin-orbit couplin g* 10 ! 11 ! f or neutral 
atoms. Several groups have also suggested the possibil- 
ity of using microwave and optical techniques in solids 
to realize topologically non-trivial effective band struc- 
tures in "trivial" materiali^MI. Topological phenom- 
ena have been identified for st rongly driverll^fl' or dis- 
sipative quantum systems^^H as we ll. These advances 
motivate the detailed study of topological phenomena in 
periodically-driven systems. 

Here we focus our attention on topological fea- 
tures of the single-particle properties of two-dimensional 
translationally-invariant tight binding systems. In the 
static case, where the Hamiltonian is constant in time, 
the topological properties of these systems are well un- 
derstood. When no additional symmetries are present, a 
complete topological characterization is provided by the 



set of values of an integer topolog ical inv ariant , the Chern 
number, evaluated for each bancplM 2 ^. 

One of the most striking applications of this topological 
characterization is to the edge physics of these systems. 
In a finite geometry with an edge, a two-dimensional sys- 
tem may support chiral edge modes which propagate at 
energies within a bulk band gap. Remarkably, the val- 
ues of the Chern numbers can be used to predict the net 
number of chiral edge modes traversing each bulk gap 
(counted according to their chirality). This edge state 
count is significant because it is a robust property which 
is independent of the fine details of the system and of the 
edge. In the case of the quantized Hall effect, the guaran- 
teed existence of these modes provides a powerful frame- 
work for explaining a variety of complex phenomena^. 

Given the power of this approach, we seek to generalize 
these results to periodically driven systems. In the driven 
case, the analogue of the energy spectrum is the Floquet 
spectrum - the set of eigenvalues of the time evolution 
operator, evaluated over one complete cycle of driving. 
The eigenvalues of the unitary evolution operator are unit 
modulus complex numbers, which are defined on a circle. 
We express them in terms of a "quasi-energy" e as e~ teT , 
where T is the driving period. The spectrum is thus 
27r/T-periodic in e. 

There are many close analogies between static and 
periodically-driven systems. In particular, the Floquet 
spectrum can be organized into quasi-energy bands with 
associated Chern numbers, just as in the static case. 
Also, driven systems may exhibit chiral edge modes when 
defined in a finite geometry with a boundary. Despite 
these similarities, however, the Chern numbers employed 
in the static case do not give a full characterization of 
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FIG. 1. Topologically protected Floquet edge modes in a two- 
band system with trivial Chern indices, C — 0. This situation 
is made possible due to the periodicity of the quasi-energy e: 
in addition to the edge modes which cross the gap near e = 0, 
a second branch of edge modes crosses through e = n/T, thus 
connecting the top of the upper Floquet band to the bottom 
of the lower Floquet band and closing the quasi-energy cycle. 
The Chern numbers correctly yield the differences between 
the numbers of chiral modes above and below every band, but 
cannot be used to uniquely determine the edge state spectrum 
of a periodically-driven system. 



the topological properties of periodically driven systems. 
In particular, for driven systems, these invariants do 
not uniquely determine the number of chiral edge modes 
within each bulk band gap. 

Recently, new types of edge modes, which cannot be 
accounted for using the invariants developed for static 
systems, were discovered in the Floquet spectra of one 
and two dimensional periodically-driven system^^22l. A 
schematic picture showing an example of such "anoma- 
lous" edge modes in a two-band two-dimensional system 
is shown in Fig. [T] Here, two sets of co-propagating chi- 
ral edge modes are found, despite the fact that the Chern 
numbers associated with both bands are zero. This sit- 
uation cannot occur in static systems, where such edge 
modes are only produced at the boundaries of systems 
characterized by non-zero Chern numbers. 

In this paper, we develop a more complete understand- 
ing of this behavior, and identify appropriate topologi- 
cal invariants (winding numbers) for characterizing these 
new phenomena. In contrast to the Floquet-band Chern 
numbers, which only depend on the evolution operator 
evaluated over one complete driving cycle, the winding 
numbers utilize the information in the evolution oper- 
ator for all times within a single driving period. We 
show that the winding numbers fully determine the chiral 
edge mode counts in each of the Floquet gaps. Thus our 
construction provides a complete "bulk-edge correspon- 
dence" for periodically-driven, two dimensional, single 
particle systems. 

The paper is organized as follows. Sec.|T]gives a general 
discussion of the bulk edge correspondence for both static 
and periodically-driven systems. In Sec. [il] we introduce 
a simple model of a strongly driven system, similar to the 
one considered in Ref. HH1 which supports a Floquet band 
structure like the one depicted in FigJTJj. Then in Sec. Ill 



we provide a general analysis and construct the winding 
numbers of the time- dependent evolution operator, which 
fully characterize the topological features of the system. 
In Sec. |IV| we give a complementary approach for deriv- 
ing the edge state spectrum, which utilizes an analysis 
in the frequency domain. We demonstrate this approach 
for the case of weak, harmonic driving, as considered in 
Refs. [T^HIH Finally, in Sec. [V] we discuss possible ex- 
perimental realizations of these new topological spectra, 
and prospects for future investigation. Technical aspects 
of the derivations are provided in the appendix. 



I. BULK-EDGE CORRESPONDENCE IN 
STATIC AND DRIVEN SYSTEMS 

One of the main goals of this paper is to explain 
how the standard correspondence between the edge state 
spectrum and the values of bulk topological invariants 
is modified for periodically driven systems. Toward this 
aim, we begin this section with a brief review of Flo- 
quet band theory in order to define what we mean by the 
"band structure" of a periodically driven system. We 
then compare and contrast the static and driven cases, 
and give a concrete example that demonstrates the situ- 
ation where intuition from the static case breaks down. 

In general, the evolution of a system governed by a 
time-dependent Hamiltonian H(t) may be quite compli- 
cated. However, in the case where the Hamiltonian de- 
pends periodically on time, H(t + T) = H(t) for some 
driving period T, Floquet theory provides a powerful 
framework for analysis. In analogy with the usual expan- 
sion in terms of stationary eigenstates of a static Hamil- 
tonian, here the evolution is conveniently described in 
terms of a basis of Floquet states. These states are solu- 
tions to the time-dependent Schrodinger equation of the 
form \tp(t)) = |$(t))e- iet , where + T)) = |$(t)}. 
Under the action of the evolution operator U(T) over 
one complete period of driving, each Floquet state is 
mapped onto itself up to a phase: \ij){T)) = U(T)\%f;(0)} = 
e~ leT \ij}({))) . In a stroboscopic sense, the Floquet states 
play the role of stationary states for U(T). The param- 
eter £, called the quasi-energy, is uniquely defined up to 
integer multiples of lj = 2ir/T: any given solution with 
quasi-energy e can also be associated with a quasi-energy 
e = e + pu), where p is an integer, through the relation 
|i>(i)) = e lpult \®(t)). Similar to the crystal momentum of 
a system with discrete translational symmetry, the quasi- 
energy can be thought of as a periodic variable defined on 
a quasi-energy Brillouin zone —ir/T < e < ir/T. When 
a system has both discrete time and spatial translation 
symmetries, the Floquet states are labeled by e and the 
crystal momentum k. The Floquet spectrum then con- 
sists of a set of bands {e„(k)}, where n is a band index. 

We now explore the ways in which the topological prop- 
erties of Floquet bands are similar to, and different from, 
those of the conventional bands of static systems. Con- 
sider first a static (non-driven) two-band system, such as 
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the Haldane modeP^, which exhibits a well-defined band 
gap. A Chern number can be defined for each band n by 
integrating the Berry curvature associated with the two 
component spinor eigenstates {|u„(k))} over the Bril- 
louin zone^. If the parameters are such that the Chern 
numbers of the two bands are non-zero, we will find chiral 
edge modes traversing the band gap when the system is 
defined in a finite geometry with an edge. The existence 
of these edge modes is "protected," in the sense that the 
modes cannot be destroyed by any continuous changes of 
parameters unless the band gap closes. 

Mathematically, the Chern number of a band is equal 
to the difference between the numbers of chiral edge 
modes entering the band from below and exiting above. 
Because the spectrum of a static system is bounded, there 
can never be any chiral edge modes extending beyond the 
bottom of the lowest band, or above the top of the upper- 
most band. Therefore, in a static two-band system, the 
Chern number uniquely defines the net number of chiral 
edge states crossing the band gap. 

Now consider a periodically driven system. We can 
define a Chern number for each Floquet band n in the 
same way as in the static case: the Chern number is 
obtained by integrating the Berry curvature associated 
with the Floquet states {\ip n (\c,t})}, evaluated at any 
fixed time t, over the crystal momentum Brillouin zone 
(see Sec. IIIC for a precise definition). For concreteness 
we typically evaluate the Chern numbers using the Flo- 
quet state wave functions at time t = 0. However, the 
Chern number is in fact independent of the time at which 
the Floquet states are evaluate d 12 1 18 1 . 

Like the static case, the Chern number of a Floquet 
band is equal to the difference between the numbers of 
chiral modes above and below a given Floquet band. 
However, due to the periodicity of quasi-energy, the Flo- 
quet spectrum is not bounded. In particular, in a two- 
band system, a chiral edge mode that extends out above 
the top of the "upper" band can pass through the quasi- 
energy Brillouin zone edge at e = tt/T, and enter the 
bottom of the "lower" band, from below. Using the rule 
that the Chern number gives the difference between the 
numbers of edge states above and below a band, we see 
that chiral edge states can be found between the two 
bands of a system with zero Chern numbers, provided 
that a "winding" edge state also connects them through 
the quasi-energy zone edge. 

To understand how such a situation can arise in a 
periodically-driven system, it is useful to consider the 
following thought experiment. Suppose we start with a 
system in which the Floquet bands have vanishing Chern 
numbers, and in which there are no chiral edge modes in 
a finite geometry. The band structure of a strip, as a 
function of the momentum fey along the strip, will re- 
semble the one shown schematically in Fig. [3k. Note 
that, as discussed above, there are two gaps to consider 
between the two bands: the gap centered at zero quasi- 
energy and the gap centered at e = tt/T. Now imagine 
that we tune some parameters in the Hamiltonian. As 



parameters are changed, one of these gaps, e.g. the one 
at e = tt/T, may close and reopen, in such a way that 
the Chern numbers of both bands become non-zero. Af- 
ter the reopening of the gap, the quasi-energy spectrum 
will look like the one in Fig. |3]d, with chiral edge modes 
crossing the e = tt/T gap. As parameters are varied fur- 
ther, the e = gap may close and reopen, bringing the 
Chern numbers of both bands back to zero. However, 
the chiral edge modes in the e — tt/T gap cannot disap- 
pear during this process, since the e = tt/T gap remains 
open throughout it. Therefore, after reopening the e = 
gap, another chiral edge mode must appear around e = 
(Fig. [3J;) such that the difference of the number of chiral 
edge modes below and above each band is zero. Thus we 
may obtain the situation where chiral edge modes exist 
despite all Chern numbers being zero. 



II. MODEL OF ANOMALOUS EDGE STATES 

In this section we illustrate the breakdown of the tra- 
ditional bulk-edge correspondence for a Floquet system 
by analyzing an explicit two-dimensional tight-binding 
model. Consider a tight-binding model on a bipartite 
square lattice, with hopping amplitudes varied in a spa- 
tially homogeneous but time-periodic way as shown in 
Fig(2k. In addition to the cyclic modulation of hop- 
ping amplitudes, a constant sublattice potential 5ab = 
£a—£b which distinguishes between A and B sites (filled 
and empty circles, respectively) may be applied. The 
system's evolution is generated by the time-dependent 
Hamiltonian 

^) = E( c L c U)^(k,t)(^) a) 

k \ > / 

4 

ff(M) = -E J «W (e lb "'^ + +e- ib »- k <7-) 



Here creates a particle in a Bloch state with crystal 
momentum k on sublattice a = {A, B}, and J n con- 
trols hopping from each B site to its neighboring A site 
along the highlighted bond in step n, shown in Figj2k. 
The Pauli matrices a z and g ± = (a x ± ia y )/2 act in 
the sublattice space. The vectors {b^} are given by 
b x = b 3 = (a, 0) and b 2 = — b 4 = (0, a). 

One driving cycle consists of five equal length seg- 
ments, of duration T/5, where T is the driving period. 
During the n-th segment of the cycle, where n = 1, . . . , 4, 
the hopping amplitude J n is set to a value J, while the 
other three hopping amplitudes are set to 0. During step 
5 all hopping amplitudes are set to 0, but the sublattice 
potential is still allowed to act. This "holding period" is 
needed in order to ensure that the model has a sufficiently 
rich phase diagram to fully illustrate the topological clas- 
sification that we will develop below. The driving proto- 
col is inherently chiral, as the cycle can be executed in 
two inequivalent patterns (1 — 2 — 3—- ■ • or 5 — 4—3—- • ■ ). 
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FIG. 2. Two-dimensional tight-binding model exhibiting 
anomalous Floquet edge modes, a) Driving protocol. Hop- 
ping amplitudes are varied in a spatially-homogeneous but 
chiral, time-periodic way. A sublattice potential Sab differ- 
entiating the two types of sites in the bipartite lattice (filled 
and open circles) is applied throughout. During four of the 
five equal-length phases of the cycle, hopping along one of the 
four distinct bond types is allowed, with amplitude J (bold 
lines). During the fifth phase, all hopping amplitudes are zero, 
while the sublattice potential remains constant, b) In the sim- 
ple case JT = 5ir/2, Sab = 0, particles in the bulk move in 
closed trajectories encircling plaquettes, and the Floquet op- 
erator is the identity. Chiral edge modes propagate along the 
boundaries, c) Floquet spectrum for the case described in b). 



For an infinite system (or a finite system with periodic 
boundary conditions), translational invariance can be ex- 
ploited to reduce the problem of finding the Floquet op- 
erator for this system to multiplying a small number of 
Pauli matrices. The simplest case, illustrated in Figs[2]3 
and c, occurs when JT/5 — tt/2 and Sab — 0. Here, 
a particle moves with probability 1 between neighboring 
sites during each hopping step of the cycle. As shown 
by the light blue trajectory in Figj2]3, over one com- 
plete driving cycle each particle makes a loop around 
a plaquette and returns to its initial position. There- 
fore the bulk Floquet operator for this case is simply the 
identity. The Floquet spectrum features two degenerate 
bands collapsed at quasi-energy zero. Because the bands 
are fully collapsed, none of the standard invariants for 
two-dimensional systems can take non-trivial values. 

What happens in a finite system with an edge? 
Naively, it appears that the Floquet operator U(T) = 1 
describes the trivial stroboscopic dynamics of a system 
with an effective Hamiltonian H e g = 0. However, us- 
ing the strip geometry shown in Fig[2j} it is straightfor- 
ward to check that the system supports chiral propagat- 
ing edge modes localized on the boundaries. Over each 
complete driving cycle, a particle moves by one unit cell 
along the edge. These modes appear in the spectrum 
as two linearly-dispersing branches with group velocities 
de/dk = ±1/T, see Figj2{;. 

Now consider the time-reversed cycle (5 — 4 — 3 — • • • ) , 
which has the opposite chirality. Here, although par- 
ticles circle around the plaquettes in the anti-clockwise 



direction, the Floquet operator (which is still identity) 
is the same as that of the original cycle. Thus the in- 
formation about the circulation direction, which is con- 
tained in the bulk evolution operator for intermediate 
times within the driving cycle, is absent in the Floquet 
operator. Physically, this information is of crucial im- 
portance, however, as it determines the propagation di- 
rection of the edge states when the system is terminated. 
Thus it is clear that a full topological characterization of 
these new phases requires an invariant defined in terms 
of the evolution operator evaluated throughout the entire 
driving period. 



III. CONSTRUCTION OF THE INVARIANT 

In this section we define an integer invariant which can 
be used to correctly predict the edge state spectrum for 
a two dimensional periodically driven system. We begin 
by setting up the problem. Let us consider a periodi- 
cally driven, tight binding model defined on a 2D lattice 
with N sites per unit cell - i.e. a generalization of ([I]). 
Suppose that the hopping amplitudes are translationally 
invariant and have finite range. Working in k-space, the 
Hamiltonian can be written as 

= c k,a^«'(k,t)c k , Q /, (2) 

kaa' 

where a, a' — 1,...,N label the sites in each unit cell and 
k lies in the first Brillouin zone. All of the bulk properties 
of the system are encoded in the N x N Hermitian matrix 
H(h,t). In particular, the bulk time evolution operator 
can be computed as 

C/(k,i) = Texp ^-ijf dt ' H(k,t')J , (3) 

while the bulk Floquet operator corresponds to the spe- 
cial case U(k, T). 

Let us suppose that the Floquet spectrum has a gap 
extending over some finite interval [e — Ae, e+Ae]. Then, 
if we define the model in a geometry with an edge, any 
Floquet eigenstates with eigenvalues lying in this interval 
must be localized near the boundary. These eigenstates 
correspond to edge modes. In analogy with the time- 
independent case, we expect that the number of these 
edge modes - counted with a sign corresponding to their 
chirality - is completely determined by the bulk time evo- 
lution operator U. Therefore, one should be able to com- 
pute the number of chiral edge modes at quasi-energy e 
given only ?7(k, t). We will now construct an explicit 
formula, defined in terms of {U(k,t),e}, which gives ex- 
actly this number. (Here {U(k, t)} denotes the set of 
evolution operators for all times within a driving period, 
< t < T). 
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A. The case of a trivial Floquet operator 

Our construction proceeds in two steps. First, we con- 
sider the special case where the (bulk) Floquet operator 
is simply the identity, i.e. U(k, T) = 1 for all k, as in 
the example described in Sec. [TTJ We then generalize to 
arbitrary Floquet operators. 

If the Floquet operator is the identity, then the Floquet 
spectrum is gapped everywhere except at e = 0. There- 
fore, to each {U(k, t)} we should be able to unambigu- 
ously associate an integer n e dge which counts the number 
of edge modes propagating across the gap (i.e. winding 
around the quasi-energy Brillouin zone). Furthermore, 
this integer must be invariant under smooth deforma- 
tions of U that preserve the condition U(k, T) = 1, since 
the number of edge modes cannot change under a defor- 
mation in which the gap remains open. 

Purely mathematical considerations suggest a natural 
guess: notice that U is periodic in k x ,k y and t (since 
U(0) = U(T) = 1 by assumption). Thus, U defines 
a map from S 1 x S 1 x S 1 — > U(N). Such maps are 
known to be classified by an integer topological invariant 
or "winding number" defined by^ 



W[U] 



1 

8^2 



■ Tr (U^dtU ■ [U- 1 d km U,U- 1 d kll U\) . (4) 

It is natural to guess that n c d gc is related to W[U] in some 
way. In appendix [A] we show that this guess is correct - 
in fact, the two integers are identical: 



n edge = W[U}. 



B. The general case 



(5) 



Next we consider the general case, where the Floquet 
operator £7(k, T) can be arbitrary. We would like to com- 
pute the number of edge modes at a quasi-energy value 
e lying within a Floquet gap. One way to do this is to 
reduce the problem to the previous case. The idea is to 
construct another time evolution operator U s satisfying 
several properties. The first property is that U e has a 
trivial Floquet operator: ?7 e (k, T) = 1 for all k. The sec- 
ond property is that there exists a one-parameter family 
of evolution operators {U s : s £ [0,1]} that smoothly 
interpolates between U and U e : 

E/ a= o(M) = tf(k,t) , U s=1 (k 1 t) = U e (k,t). (6) 

Finally, and most importantly, we require that the in- 
terpolation U s (k, T) maintains a gap around some quasi- 
energy value e s which changes smoothly in s and satisfies 
£ s =o = e an d £«=i = tt/T. 

If we can construct an evolution operator U e with these 
properties then we can immediately compute the number 
of chiral edge modes of U at quasi-energy e: 



The validity of this identification comes from the fact that 
the gap doesn't close during the interpolation process (by 
assumption). Therefore, the number of edge modes of U 
at quasi-energy e must be the same as the number of edge 
modes of U e at quasi-energy tt/T. The latter quantity is 
then given by W[t/ e ], using the formula in Eq. ([5]). 

All that remains is to construct an appropriate U e and 
a corresponding interpolation. There is some arbitrari- 
ness here, since U e is far from unique. The result, how- 
ever, will not depend on our particular choice. We will 
define U e by 



U e QL,t) 



C7(k, 2t) if < t < T/2 

V E (k, 2T - 2b) if T/2 <t<T, 



(8) 



where 



14(k, t) = e-^^ , H cS (k) = I log [/(k, T). (9) 

Here, we choose the branch cut of the logarithm to lie 
along the direction e~ ieT . That is, we choose a branch 
with 



loge- i£T+40 ~ =-i-eT 
\oge' ieT+M+ =-%-eT -2-kL 



(10) 



This choice of branch is important, and provides the only 
dependence of U e on e. 

Physically, V e can be viewed as a trivial "return map" 
which is used to connect the Floquet operator U(k,T) 
to the identity. The quantity H e ff appearing in the ex- 
ponent plays the role of a static effective Hamiltonian 
which generates a bulk time evolution that, when exam- 
ined stroboscopically at integer multiples of the driving 
period T, is identical to that of f/(k, T}p3. By concate- 
nating the driving cycle with an evolution generated by 
—H e g, the net effect of the combined evolution becomes 
trivial. 

We now check that the above definition of U e satisfies 
all of our requirements. It is clear that U E (k, T) = 1; the 
crucial question is to find an appropriate interpolation 
connecting U and U £ . The following interpolation does 
the job: 



U s (k,t) 



C/(k, (1 + s)t) 



if 0<t<T/(l + s) 



V e {k,2T-(l + s)b) if T/(l + s) < t < T. 



Indeed, using Eq. ^ it is easy to check that the Floquet 
gap remains open around e s = (1 — s)(s + n/T) — tt/T, 
using the fact that 



U s (k,T) = U(k,T) 1 



(11) 



«<edge(e) = W[U £ 



(7) 



Also, we can see that e s= o = £ & n d £ s =i = tt/T (mod 
2tt/T). Thus, the interpolation U s satisfies all of the 
conditions listed above. 
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Relation with Chern number 



As discussed in Sec. |TJ the Chern number of a Floquet 
band n is defined by integrating the Berry curvature of 
the Floquet eigenstates |^„(k, f)) over the crystal mo- 
mentum Brillouin zone, at a fixed time t: 



r 



1 

27T 



dk x dk y (V x A n ), 



(12) 



where A n — (VVi(k, t)\i^\tp n (k, t)). Equivalently, C can 
be written as 

Cn = ^riJ dkxdky Tr ( Pn ' ^ P ^ dk v P ^) ' ( 13 ) 

where P„(k) = \ip n (k,t)){tp n (k,t)\ is a projector onto the 
Floquet eigenstate \ip n (k,t)). 

Not surprisingly, there is a close mathematical relation- 
ship between the winding number W[i7 e ] defined above 
and the Chern numbers {C n } of the Floquet bands. This 
relationship becomes clear when one considers the differ- 
ence between winding numbers evaluated at two different 
quasi-energies. More specifically, it is possible to show 
that 



W[U e ,] - W[U e ] = c e£/ 



(14) 



where C ee i denotes the sum of the Chern numbers of all 
Floquet bands that lie in between e and e' (see appendix 
[b] for a derivation) . This identity is very natural from a 
physical point of view. Indeed, identifying VF[C/ e ] with 
^edge(e), Eq. (14) is simply the statement that difference 



between the numbers of edge modes at two quasi-energies 
e and e' is equal to the total Chern number of the inter- 
mediate bands. 



D. Phases of the modulated square lattice model 

We now briefly illustrate the utility of the winding 
number by examining the phase diagram of the model in- 
troduced in Sec.[n] Away from the special point Sab = 0, 
JT/5 = 7r/2, the Floquet operator is not equal to 
the identity, and the return-map construction described 
above must be used to compute the winding numbers 
for the Floquet gaps around quasi-energy values and 
7r. Using direct numerical calculation, we solved for Flo- 
quet spectra in both periodic and strip geometries (see 
e.g. Figs. [3^,-c), and mapped out the approximate phase 
diagram, shown in Fig. |3ji. Example spectra from each 
of the three phases confirm the correspondence between 
the winding number and the edge state count in each 
gap. Figure |3j; provides an explicit example of the phe- 
nomenon discovered in Ref.[18]and shown in Fig.[l} chiral 
edge modes appear for the finite system, despite the fact 
that the Chern numbers of both bands are zero. 
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FIG. 3. Phases of the square lattice model in Sec. [IT] a-c) Ex- 
ample Floquet spectra for each of the three phases, calculated 
in a strip geometry (unit cell 2a), with sublattice potential 
Sab = 0.5n/T and hopping amplitudes a) J = 0.5 n/T, b) 
J = 1.57I-/T and c) J = 2.5tt/T. The winding numbers W 
and Wtt correctly yield the numbers of edge states in the two 
gaps. The asymmetry in ka is due to the breaking of inversion 
symmetry by the sublattice potential Sab- d) Phase diagram 
indicating the winding numbers Wq and W n , calculated in the 
gaps at quasi-energies and n/T, and Chern number C of the 
upper band. Note the existence of two topologically-distinct 
phases with C = 0. 



IV. FREQUENCY DOMAIN FORMULATION 

In this section we discuss an alternative approach for 
deriving the edge state spectrum of the Floquet opera- 
tor, and its relation to the Chern numbers of the bulk 
Floquet bands. This approach is based on an analysis in 
the frequency domain instead of the time domain. 



A. Repeated zone analysis 

We start from the Schrodinger equation with Hamil- 
tonian ^ for the Floquet state in band n, with crystal 
momentum k. Using a basis of states labeled by a, we 
write \ip n (k,t)) = X)„=i Vw«(k, ^cj^jO), where |0) is the 
vacuum. The amplitudes {ip na (k, t)} evolve according to: 

JV 

id t 1p na {k,t) = ^2 H a a >(Kt)Tpna'{Kt)- (15) 



Employing the Floquet theorem, we write 



i>na(Kt) 



_ p -ie„(k)i 



(16) 
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where u> — 2tt/T. Below we suppress all k indices for 
notational simplicity. The coefficients tpf^ satisfy the 
(time-independent) eigenvalue equation 



a' ,m' 



H 



mm , J m ') 
act' t na f 



(17) 



where the "Floquet Hamiltonian" H™™ is given by 



HI 



moj6 aa ,6 mm , + - / dte-^-^H^it). 
1 Jo 

(18) 

For each k, Eq. (17) has solutions throughout — oo < 
e„ < oo. However, as discussed in Sec. [TJ if e„ is an 
eigenvalue of Eq. ( [l7| corresponding to the eigenstate 

with amplitudes {ipna}, then e n = e„ + pui (where p 
is any integer) is also an eigenvalue corresponding to an 

In 



eigenstate with amplitudes given by 
fact, as seen by direct substitution into Eq. (|T6|), all of 



(m) 
^Pna • 



these solutions correspond to the same time-dependent 
solution of the Schrodinger equation. Therefore, the Flo- 
quet states are uniquely and completely parametrized by 
quasi-energies in the "first quasi-energy Brillouin zone," 
— 7r/T < e„ < w/T. Equation (17) is the temporal ana- 
logue of a "repeated zone" scheme for conventional band 
structure calculations. 

To illustrate the structure of Eqs. (17) and (18), we 
now consider the case where 



H(t) = H Q 



Ae iut 



Ate 



(19) 



In this case, the matrix H™™ has the block tri-diagonal 
form shown schematically in Fig. [4^,, where each block is 
an N x N matrix. As noted above, it suffices to solve for 
the Floquet bands within the first quasi-energy Brillouin 
zone, —tt/T < s n < tt/T. In the limit of weak driving 
(small A), this can be done in perturbation theory. 

As sketched in Fig. |4Jd, the zeroth-order spectrum of 
the Floquet Hamiltonian ( 18 ) consists of an array of 



copies of the original spectrum of Ho, shifted up and 
down by integer multiples of the drive frequency lo. In 
the Floquet picture, the harmonic driving (Ae lujt + h.c.) 
induces hopping between levels with m-values differing 
by one. Processes where m increases (decreases) by one 
correspond to transitions accompanied by the absorption 
(emission) of a photon from the driving field. When the 
range A of the spectrum of Hq is greater than w, resonant 
transitions (associated with degeneracies between levels 
with neighboring values of m) may occur for some values 
of k. To leading order, degenerate perturbation theory 
shows that the driving (Ae lut +h.c.) opens gaps at these 
resonances (see FigM 



B. 



Truncation of the Floquet Hamiltonian and 
edge states 



An important property of Floquet eigenstates is that 
they are localized in m, decaying rapidly for \m — mo\uj 3> 



a) 



U(k) 




b) 



A 



LO 



-2 -1 





m 



FIG. 4. Floquet Hamiltonian and level diagram in the re- 
peated zone scheme, a) When the driving A(t) consists of a 
single harmonic with a ngu lar frequency w, the Floquet Hamil- 
tonian H(k), see Eq. (18 1, assumes a block-tridiagonal form. 
Each block is an iV x N matrix, where iV is the number of 
bands in the unperturbed Hamiltonian Hq. The off-diagonal 
blocks labeled by A and A* describe transitions accompa- 
nied by the absorption and emission of a driving field photon, 
respectively, b) Schematic Floquet level diagram. The iV- 
level spectrum of Ho(\c), which spans an energy range A, is 
copied and rigidly shifted by moj for each value of m. The 
harmonic driving A induces transitions between neighboring 
copies with m-values differing by 1. Analogous to Wannier- 
Stark states, the Floquet eigenstates in the first quasi-energy 
Brillouin zone (gray shaded region) are localized in m, with 
amplitudes highly suppressed for \m\ S> A/w. 



A, where mo is the center of a given localized state. This 
localization in frequency space is analogous to the well- 
known localization of Wannier-Stark states in real space. 
Starting from this observation, we now describe an alter- 
native method for computing edge state spectra of Flo- 
quet systems - complementary to the winding number 
approach of Sec. |III| 

The first step is to truncate the Floquet Hamiltonian 



(18) so that m and m! run over a large but finite range, 
— M < m, m' < M, where M is much greater than the 
frequency-space localization range of the Floquet states. 
Making use of the localization of the Floquet states in 
m, we note that the spectrum of the truncated Floquet 
Hamiltonian will be a good approximation to the exact 
result within the first few quasi-energy zones centered 
around m = 0. In particular, if we consider a geome- 
try with a boundary, then the edge state spectrum of the 
truncated Hamiltonian will closely match the exact Flo- 
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quet edge state spectrum within these central few quasi- 
energy zones. 

The key observation is that there is actually a straight- 
forward way to determine the edge state spectrum us- 
ing Chern numbers of the bands of the truncated Flo- 
quet Hamiltonian. Interpreting the truncated Floquet 
Hamiltonian as the static Hamiltonian of some new 
(2M + 1) x TV-band system, we apply the standard bulk- 
edge correspondence to deduce that the net number of 
chiral edge modes (counterclockwise minus clockwise) 
crossing any particular gap is given by the sum of the 
Chern numbers of all bands below this gap. In this way, 
we have a simple procedure for deriving the edge state 
spectra of Floquet systems: to find the number of chiral 
edge modes crossing a given quasi-energy gap within the 
first quasi-energy Brillouin zone, i.e. within the interval 
—tt/T < e < tt/T, we truncate the Floquet Hamiltonian 
and then count Chern numbers of the resulting bands 
below the gap in which we are interested. Note that, 
as in Sec. 



IIIC 



the difference between the numbers of 
chiral edge modes at quasi-energies e and e' within the 
first quasi-energy zone is equal to the sum of the Chern 
numbers of all Floquet bands between s and e' . 

As an example, we study the case of harmonic driving 
introduced in Eq. ( 19 1. We take H to be a general two- 
band Hamiltonian, 



tf (k) = e(k) + d(k)-<T, 



(20) 



where cr = (a x , a v , a z ) is the vector of Pauli matrices, 
and d(k) is a three dimensional vector. The bands are 
characterized by Chern numbers C — ±Co, as indicated 
in Fig. [5^,. For simplicity we assume that the periodic 
drive creates a single resonance between the valence and 
conduction bands of Hq. This resonance occurs simulta- 
neously for all k-values which satisfy uj — 2|d(k)|. We 
further assume that this resonance condition is satisfied 
along a single closed curve in the crystal momentum Bril- 
louin zone (see Fig. [6^i). 

Consider the co rresponding Floquet Hamiltonian 
, see Eq. (18), which is truncated to the range 
— 1 < m,m! < 1. The band structure of the truncated 
Floquet Hamiltonian is shown schematically in Fig. [5] 
For A = 0, the valence band of block m and the con- 
duction band from block m — 1 are degenerate on the 
resonance curve (appearing as two points on the sec- 
tions shown in Fig. [5]) . The driving terms mix the two 
bands, opening gaps at the crossing points. The resulting 
bands have Chern numbers Cf and —Cf, as indicated in 
Fig. [5Jd. Due to the truncation, however, two isolated 
bands are left, at the top and bottom of the spectrum. 
These bands do not participate in any resonances, and 
their Chern numbers remain equal to ±Co, see Fig.JS}). 

Examining the gap centered at quasi-energy —tt/T, we 
see that the Chern numbers of all bands below it sum to 
Co + Cf- Therefore we expect to find Co + Cf chiral 
edge modes spanning this gap. For the same reason, the 
gaps centered at quasi-energy values e = and e = tt/T 
support Cq and Co + Cf chiral edge modes, respectively. 



y^~m il^V...... 



-C = -Co 




b) 



- 



c 



-C F 





FIG. 5. Floquet bands for a two-band model in the trun- 
cated repeated zone scheme. Representative one-dimensional 
cuts through the two-dimensional crystal momentum Bril- 
louin zone are shown, a) For A = 0, the unperturbed bands 
of Hq are each copied and shifted up and down one time, 
m = —1,0, 1. The bands have Chern numbers C = ±Ct>. The 
crossings between bands with m-values differing by 1 corre- 
spond to points in the Brillouin zone where the driving field 
resonantly couples the two bands, b) For A/0, the driving 
field opens avoided crossings at the resonance points. After 
coupling, the Chern numbers of the crossed bands are changed 
to C — ±Cf- In the truncated scheme, the appearance of 
Floquet edge states at a particular quasi-energy in the first 
quasi-energy Brillouin zone (gray region) can be predicted by 
adding up the Chern numbers of all bands below it. 



The value of Cf which is obtained in any particular 
model depends on the details of Ho(k), A, and oj. In 
particular, in Sec. |IV C| below and in Fig [6] we describe 
a specific model which gives Cf = 0, while nonetheless 
exhibiting robust chiral edge modes. 

Importantly, the method outlined above for determin- 
ing the number of chiral edge modes in each gap of the 
Floquet spectrum can be applied to more general period- 
ically driven systems, and is not restricted to the weakly 
driven, single harmonic limit discussed in the above ex- 
ample. It works provided that one keeps sufficiently 
many copies of Hq in the truncation, and that the driv- 
ing power spectrum decays sufficiently rapidly for higher 
harmonics. If the driving is strong, the relationship be- 
tween the Floquet bands and the original bands may be 
more complicated than in the example considered above, 
but the results of this procedure are still guaranteed to 
converge for large enough M. 



C. Anomalous edge states in a weakly-driven 
two-band lattice system 

We consider a two-band model with a Hamiltonian of 



the form (20 1, with e(k) = and with the components 
of d(k) given by: 



d x (k.) = asin(fe x ), d v (\i) = as'm(k y ) 
d z (k) = (/i — J) — 26(2 — cos(fca;) — cos(fc J/ )) 
+Jcos(fe x ) cos(ky). 



(21) 
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In the parameter regime n, b, J > and fi/b < 4, the un- 
perturbed bands of Hq have Chern numbers Cq = ±1, as 
indicated in Fig.[6ji. This result can easily be obtained by 
direct integration of the Berry curvature using Eq. (12 1. 



a) -1 



More simply, recall that for a general two-band system 
of the form (201, the Chern number is determined by the 



degree of the map d(k) = d(k)/|d(k)| from the Brillouin 
zone to the unit (Bloch) sphere S 2 , i.e. by the number 
of times this map "covers" S 2 . The result Co = ±1 fol- 
lows simply by noting that the vector d(k) maps exactly 
one point in the Brillouin zone to a point close to the 
north pole of the unit sphere, as indicated by the colored 
shading in Fig. [6^, (see Ref . l29l for more details). 

We now add a time-dependent perturbation of the form 
A(t) = AoCTj cos(arf). The hopping J and frequency ui 
can be tuned such that only a single resonance occurs 
in the Brillouin zone, on a closed curve around k = 
as shown in Fig. [6^,. Note that here we take parame- 
ters to provide the simplest situation, containing only a 
single resonance, for the purposes of clearest illustration; 
similar considerations apply in other parameter regimes 
allowing for more resonances. 

To see that this driving results in "trivial" Floquet 
bands with Cp = 0, we perform degenerate perturbation 
theory in the off-diagonal matrix elements of the Flo- 
quet Hamiltonian, Eq. ( |17[ ). Consider an isolated pair of 
crossing bands, such as the red and blue bands labeled 
m = —1 and m — in Fig. [5^, respectively. To lowest 
order in Ao , the corresponding two new hybridized bands 
are described as the eigenstates of an effective Hamilto- 
nian 



H 



(o) _ 

eff — 



|d(k)| - w/2)d(k) • tr + A(k) • a, (22) 



Note that H$ is in 



where A(k) = A [z - <4(k)d(k) 
fact a leading-order approximation to the effective Hamil- 
tonian introduced in Eq. ([9|. As indicated in Fig.[5j}, the 

bands of H^) are inverted relative to those of Hq in the 
region around k = 0. 

How do we see that these bands have trivial Chern 
indices? Consider the map from the Brillouin zone to S 2 



defined by H^) in Eq.(|22|). A necessary requirement for 



obtaining a nonzero Chern number is that every point on 
the Bloch sphere must be reached for at least one value 
of k in the Brillouin zone. For small Ao , and for w such 
that the resonance curve winds around k = 0, there is 
a region surrounding the north pole of S 2 which is not 
accessed by H^g for any points in the Brillouin zone. 
Therefore its bands must have zero Chern numbers. 

In Fig. [6bj we plot the spectrum of the truncated Flo- 
quet Hamiltonian for the model described above, defined 
in a cylindrical geometry. This figure clearly exhibits 
chiral edge states traversing each of the gaps, including 
those associated with the isolated bands at the top and 
bottom of the truncation window. As described above, 
these bands are not strongly hybridized, and are nearly 
identical to the unperturbed bands of Hq (see Fig. fflt 




FIG. 6. Anomalous edge states in a weakly-driven two-band 
system, with Hamiltonian Ho given by Eqs. (201 and (21 1. 



The Chern numbers of the unperturbed bands are Co = ±1. 
a) Spectrum of Ho with periodic boundary conditions, show- 
ing the resonance curve (purple line) in the valence and con- 
duction bands. The color scheme indicates (a z ) for the cor- 
responding eigenstates in the two bands, b) Spectrum of 
Ho for a cylindrical geometry, with periodic boundary con- 
ditions in the x-direction and open boundary conditions in 
the y-direction. c) Spectrum of the truncated Floquet Hamil- 
tonian, Eq. (181, in the cylindrical geometry. The isolated 



bands at the top and bottom of the figure do not partici- 
pate in any resonances and are nearly identical to the original 
bands of Ho- The Floquet bands have zero Chern numbers, 
while edge states traverse both gaps, centered around e — 
and e = tt/T. The parameters used are J/[i = b/fi = 1.5, 
a/fi = 4, Ao/fJ, = 1, and A /w = 0.07. 



The analysis above provides an alternative pic- 
ture for understanding how edge states can appear 
in periodically-driven systems, even when all of the 
Floquet bands have vanishing Chern numbers. The 
fact that the bands carry zero Chern numbers only 
guarantees that the numbers of chiral edge states in the 
two gaps surrounding each Floquet band are equal. In 
the truncated Floquet Hamiltonian picture, the edge 
states owe their existence to bands near the truncation 
boundaries which carry anomalous (i.e. non-vanishing) 
Chern numbers. These bands retain the history about 
the topological properties of the original bands of Hq, as 
well as the structure of the resonances which transform 
them into the Floquet bands. Note that anomalous 
edge states can be obtained also in a model when the 
unperturbed Hamiltonian ffo(k) contains only trivial 
bands, provided that additional resonances (e.g. due to 
multi-photon processes) are allowed. 
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V. DISCUSSION 

In this paper, we discussed the correspondence between 
bulk topological invariants and the edge state spectra 
of two dimensional periodically driven systems. In both 
static and periodically-driven systems, the Chern number 
of a given band is equal to the difference between the 
number of chiral edge states above and below the band. 
Because the spectrum of a static system is bounded from 
below, knowledge of the Chern numbers of all bands up 
to a particular energy is sufficient to uniquely determine 
the number of chiral edge modes at that energy. For 
driven systems, however, the spectrum is defined on a 
circle, and hence is not bounded. Knowing the Chern 
numbers of all Floquet bands is therefore not sufficient 
to determine the edge state structure. Our main result 
in this work is the construction of a new invariant which 
fully captures the edge state spectrum for periodically 
driven systems. 

The most striking example of the difference between 
the static and driven cases occurs when a driven sys- 
tem supports robust chiral edge states, even when the 
Chern numbers of all of its Floquet bands are zero. Can 
such systems be realized experimentally? In the anal- 
ysis above we introduced two lattice models: one with 
strongly modulated hopping amplitudes, and one with a 
weak driving field applied on top of a convenient band 
structure. While the first model was constructed pri- 
marily for mathematical demonstration, a class of mod- 
els similar to the second, involving weak, uniform driving, 
should be accessible in cold atomic or solid state systems. 

A promising route for implementing models of the 
second type is to utilize Cooper's optical flux-lattice 
systems^. There, weak periodic modulations of the op- 
tical fields which create the flux lattice could be used 
to produce the resonances necessary for obtaining a 
non-trivial Floquet spectrum. A modulated flux lat- 
tice with "anomalous" edge states could be achieved 
cither by starting with topologically-nontrivial unper- 



turbed bands, as in Sec. |IVC| or even by starting with 
topologically-trivial bands and employing two (or more) 
photon resonances. 

Solid state implementations are also possible. Typi- 
cally, the frequency uj will be at least a few times smaller 
than the total width A of the two central bands of interest 
(e.g. the conduction and valence bands). In this case, the 
Floquet band structure is complicated by quasi-energy 
zone folding, which leads to additional band crossings 
corresponding to multi-photon resonances. Generically, 
for w/ A <C 1, one can find bulk states coming from a suffi- 
ciently high-order folding at any quasi-energy. Therefore, 
strictly speaking, there is no gap in the Floquet spectrum 
at any quasi-energy. If there are chiral edge states in a 
gap which corresponds to a first-order (single photon) 
resonance, they acquire a finite lifetime due to scattering 
into bulk states with the same quasi-energy. However, 
these scattering processes involve high-order resonances 
(multi-photon processes). For sufficiently weak driving, 



the lifetime of the Floquet edge states is parametrically 
smaller than the size of the gap in which they reside. 
Therefore, the chiral Floquet edge states are still well- 
defined excitations . 

Throughout this work, we have focused on the topo- 
logical characteristics of the single-particle spectra of 
periodically-driven two-dimensional systems. Notably, 
we have not discussed the filling of these states for 
fermionic (or bosonic) systems. This is a very interesting 
and challenging problem of its own, which is beyond the 
scope of this work. However, future studies on the phe- 
nomenology of periodically-driven multi-particle systems 
will need to address this issue in detail. 

The invariants presented in this paper give a complete 
topological description of the Floquet band structures of 
periodically driven, two-dimensional non-interacting sys- 
tems, if no additional symmetries are imposed. In static 
systems, the presence of additional symmetries leads 
to a rich topological characterization of non-interacting 
systematise. i n the case of periodically driven systems, 
some robust topological phenomena have been demon- 
strated to be stabilized by additional symmetries such as 
an effective time reversal symmetr y^ 12 * 14 * 16 !. A full topo- 
logical classification for driven systems is yet to be devel- 
oped, and is an interesting direction for further study. 
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Appendix A: Derivation of the relation between 
edge modes and winding number 

In this section, we prove equation ([5]): we show that the 
number of Floquet edge modes is given by the winding 
number, i.e. n c d gc = W[J7], for any system in which the 
bulk Floquet operator is the identity. Our derivation is 
inspired by an analogous derivation in Ref. 1281 

The precise statement we will prove is as follows. 
Let H be a translationally invariant, periodically driven 
Hamiltonian with a time evolution operator U. Sup- 
pose that U has a trivial Floquet operator - that is, 
(7(k, T) = 1. Now, consider a cylindrical geometry with 
periodic boundary conditions in the y-direction and open 
boundary conditions at x = 1 and x = L x (here x and 
y are integers labeling lattice sites). Consider a Hamil- 
tonian H which is identical to H in the interior of the 
cylinder, but can take any form near the boundaries as 
long as it is local and translationally invariant in the y 
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direction. Let U be the corresponding time evolution op- 
erator. What we will show is that the number of edge 
modes of E/(k, T) is given by the winding number W[U]. 

The first step is to find a formal mathematical expres- 
sion for the number of Floquet edge modes, n e d gc - To 
this end, it is convenient to describe H and U in mixed 
x, k y space. In this description, we denote the Hamilto- 
nian by H xx i(k y ,t) and the time evolution operator by 
U xx '(k y ,t). Here we drop the additional indices a, a' 
corresponding to the sites within the unit cell, so that 
H XX '(k y ,t) and U xx >(k y ,t) are actually N x N matrices 
for each x, x' , k v , t. 

The time evolution operator U xx >(k y ,t) satisfies sev- 
eral important properties. The first property is that it is 
quasi-diagonal in x: U xx > (fc„, t) — ¥ as \x — x'\ — > oo. To 
see this, note that U xx > can be thought of as a propaga- 
tor for single particle states. Furthermore, the velocity 
of these states is bounded during the course of the time 
evolution: \v\ < c for some c. It then follows that U xx i 
vanishes exponentially if x, x' are separated by a distance 
greater than c • t. 

Another important property of U is that the Floquet 
operator U xx i(k y ,T) reduces to the identity matrix S xx i 
unless x and x' are near x — 1 or x — L x (i.e., within dis- 
tance c-T of one of the boundaries). This property follows 
from the same reasoning as above: as long as x and x' 
are far from the boundaries, the amplitude to propagate 
between them must be identical to that in the transla- 
tionally invariant system (where U(T) — 1). Therefore 



U xx / ( < ky 1 T^ — U xx '{ky 1 T^ — S xx 



(Al) 



unless x and x' are near one of the boundaries. 

It is illuminating to write out the matrix U(k y , T) more 
explicitly. All together, the matrix U(k y ,T) has dimen- 
sion NL X x NL X , since the lattice index x runs from 1 
to L x , and the band index a runs from 1 to N. Let us 
list the rows and columns in order x = 1,...,L X , with 
the band index running from 1 to N for each value of 
x. Then, it follows from the above two properties that 
U{k y ,T) (approximately) takes the block diagonal form 



U(ky,T) 




(A2) 



where Ui(k y ) and U^{k y ) are unitary matrices describing 
the action of U near the two boundaries and 1 describes 
the action of U in the interior of the cylinder. To see this, 
imagine we partition the lattice sites x into three groups: 
those with x < L x /3, those with L x /3 < x < 2L x /3, 
and finally those with x > 2L x /3. In the thermody- 
namic limit, U(k y ,T) will be block diagonal with re- 
spect to these three groups - up to corrections which 
are exponentially small in L x . Furthermore, U will act 
like the identity matrix 1 on the middle group of sites, 
L x /3 < x < 2L x /3. Thus, we can define the three diago- 
nal blocks {Ui(k y ), 1, U?,(ky)} to be NL x /3 x NL x /3 ma- 



trices which describe the action of U within these three 
groups of sites. 

Using the above matrix representation for U, we can 
write down a simple expression for the number of edge 
modes, n c d gc , on the x — L x edge. In particular, we 
claim that 

nedge = -^~J dk y Tr{U{k y ,T)- x d kv U{k v ,T)-Q), 

(A3) 

(A4) 
(A5) 



where Q is an operator defined by 

Qxx' — Q\P^)^ X x' i 

and g is any function satisfying 



9{x) = 



if x < L x /3 
if x > 2L x /3. 



The trace in Eq. (A3 1 is taken over the band indices a, a' 
and the site indices x.x'. 



To understand where ( A3 ) comes from, imagine we 



replace Q by the identity operator 1. Then the above in- 
tegral reduces to J dk y Tr(U~ 1 dk y U), which counts 
the total number of chiral modes propagating in the y 
direction - both at x = L x and x = 1^. Our claim 
is that, when we include the operator Q in the expres- 
sion, we effectively count the edge modes near the x = L x 
boundary, and throw out the contribution from the x = 1 
boundary. 

To see this explicitly, let us write out the matrix for Q 



using the same notation as equation ( A2 1 : 



'0 (V 
Q 2 
.0 1, 



(A6) 



Here Q2 describes the action of Q on the sites with 
L x /3 < x < 2L x /3 while and 1 describe the action 
of Q on the sites with x < L x /?> and x > 2L x /3, respec- 
tively. From the matrix representations for Q and U , we 
can see that the only contribution to Tr([/ -1 <9fc U ■ Q) 
comes from the sites near the x = L x boundary: 

Ti{U- x d ky U ■ Q) = TT(U^\k y )d ky U 3 (k y )). (A7) 

Clearly, when we integrate this expression over k y and 
take the trace, the result is simply the number of chi- 
ral edge modes localized near the x = L x boundary, as 
claimed above. 



The next step is to write (A3) as an integral over t: 

dtdk y d t Tr(U- x d ky U -Q). (A8) 



1 

~2tH 



We then add the total derivative d k Tr(— U 1 dtU ■ Q) to 
the integrand, thereby deriving 



. Qii/dkii 
2ni 1 y 
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Next we note that the integrand only receives con- 
tributions from the interior of the cylinder, since Q is 
constant near the boundaries and hence the commutator 
[Q, U~ 1 dtU] vanishes there. In the interior of the cylin- 
der, however, U is identical to U. We are therefore free 
to replace U —> U: 



/ dtdk v 

2m 1 v 



■TviU^dkU-iQ^U^dtU]). 



(A9) 



To proceed further, let A = U^d^U, B = U~ l d t U. 
Exploiting the translational invariance of A, B, we have 

Tr(A ■ [Q, B]) = A ^B x , x (g{x') - g{x)) 

xx' 

= ^2 A OsB s0 ■ (g(x + s) - g(x)). 

xs 

Next, we use the identity J2x(d( x + s ) — d( x )) — s to 
write this as 

Tr(A-[Q,B])=J2A 0s B s0 -s 



2n 



dk x A(k x )d kn: B{k x ), 



where the last equality comes from taking the Fourier 
transform. Substituting this expression into ( A9 ) , we 
have 



1 

4^ 



dtdk x dk y TriU^d^U- dk^U^dtU)). 



The final step is to massage this expression into the de- 
sired form ([5]) by adding a derivative with respect to 
k Xl k y , and t to the integrand. In particular, we add: 

- \d hm Tt{U- 1 d ky U > U^dtU) + ^d iv Tr(U- 1 d km 8 t U) 



1 



dtTriU-'dk^U). 



Adding this term and simplifying gives 



^edge — g 2 / dtdkxdky 



Tr(U- x d kx U- U^dkXJ ■ U^dtU) - (k x o k y ) 
= W[U], (A10) 



as claimed. 



Appendix B: Derivation of ( ]14| ) 



In this section, we derive relation ( 14 ) between the 



winding number W[J7 e ] and the Chern number of the 
Floquet bands. Let U(k,T) be Floquet operator with 
gaps at e,e'. We wish to compute the difference in the 
winding numbers W[f/ e '] — VT[J7 e ]. A simple way to do 



this is to note that FF[J7 e '] — W[E/ e ] is equal to the winding 
number of a map U which "glues" V £ > and V £ together 
with opposite orientations. To be precise, 



W[U e ,]-W[U e } = W[U], 



(Bl) 



where 



^ ' 1 V e >{k,2T -It) ifT/2<i<T. v ; 

To compute W[E7], we use the fact that it is invariant 
under continuous deformations of U . We note that the 
above map U can be continuously deformed into U where 



u(k,t) = v e (k,t)-v e ,(k,ty 



(B3) 



To see this, note that the following interpolation does the 
job: 



U s (k,t) 



V e (k,2t)-V e ,(k,2ty 
V £ ,{k, 2T -2t) 1 - s 



if < t < T/2 
if T/2<t< T, 



for < s < 1. 

The problem reduces to computing the winding num- 
ber W[E/]. Note that V £ and V e i differ only in the po- 
sition of the branch cut of the logarithm in definition 
([9]), which affects only the quasi-energy eigenvalues but 
not the eigenvectors. These two operators can therefore 
be simultaneously diagonalized, and commute. From the 
definition we see that 



U(k,t) = V e {k,t) -V e ,{k,t)- L = 



1 = „^p„,(k) 



(B4) 



where P ££ '(k) is the projector onto the Floquet eigen- 
states with crystal momentum k, and with eigenvalues 
between e and e' . Therefore, according to the calcula- 
tion of Appendix [Cj the corresponding winding number 
is W[U] — C ££ ', where C ££ > is the total Chern number 
of all the Floquet bands with eigenvalues between e,e'. 
Putting this all together, we deduce: 



W[U e ,] - W[U e ] = W[U] = C £ 
as we set out to show. 



(B5) 



Appendix C: Winding number for time independent, 
fiat band Hamiltonians 

In this section, we compute the value of W[U] for time 
independent, flat band Hamiltonians of the form 



H(k) 



2tt 

Y 



P(k), 



(Cl) 



where P(k) is an N x N Hermitian matrix, all of whose 
eigenvalues are either or 1 (i.e. P is a projection op- 
erator). Note that the prefactor ^ guarantees that the 
corresponding Floquet operator is trivial, i.e. U (T) = 



13 



e~ lHT = 1. For concreteness we assume that P(k) has 
M eigenvalues equal to 1 , and N — M eigenvalues equal 
to where M < N. Thus, the Hamiltonian H has M 
bands with energy —2it/T and N — M bands with energy 
0. We will show that W[U] is given by the total Chern 
number C of the bands with energy —2tt/T. 

First, we note that the time evolution operator U = 
e -|f-P(k) sa ti s fi eg the identities 



U 

U- 1 = (e 



l)P+ 1 
-1)P + 1, 



(C2) 



from which it follows that 

U-^U = ( e - 2 ^/ T - l)(e 2 * H ' T - l)Pd kx P 
+ (e 2 "*/ T - l)d kx P 

'2irt 



Here, 



a{6) = 2 - 2cos(0) , 6(0) = e if> - 1. 



(C4) 



Similarly, we have 



U- x d k U = a ( ^ ] • Pd k „P + b [ ^ | • d k ..P 



2nt 



2nt 



Combining these results, we find 



(C5) 



Tr ( [U- 1 dt U] [U- 1 d k=c U] [U~ 1 d ky U] ) = 

/ 27T7 

Tv( — P[a ■ Pd K P + b ■ d kx P}[a ■ Pd ky P + b ■ d ky P] 



Next, we make use of the identity 

(Pd.P) ■ P = (d t (P 2 ) - d t PP) ■ P 
= d t P P-diP P 
= 

to reduce the expression to 

TT([U- 1 d t U][U- 1 d k U}[U- 1 d k U}) = 

— (a + b)b-Tv(Pd kx Pd ky P). 
Substituting in the above expressions for a, b, we find 



^{[u^dtUW^dkUWu^dkU]) = 

2n? 

— (2cos(27rf/T) - 2) ■ Tr(Pd k „Pd ky P). (C6) 

We are now ready to compute the winding number cor- 
responding to U: 



d kx P. (C3) W W] 



8tt 2 J x y 
■ Tr ([U^dtUW^dkMP^d^U]) - (k x o k y ) 

1 f , 2iri 

~ 8^ _ 

■Tr(P.[d K P,d ky P\) 
1 



„ — (2cos(27ri/T)-2) 



2?ri 



dk x dk y Tr(P- [d ksc P,d ky P]) 



(C7) 



The latter expression is precisely the Chern number C, 
Eq. (13), of the bands with quasi-energy e = —2ir/T. 

It is also possible to derive this identity from simple 
physical considerations. To this end, let us calculate the 
number of Floquet edge modes n e dgo corresponding to 
U in two different ways. On one hand, we know that 
"edge = W[J7]. On the other hand, it follows from the 
special form of U that the number of Floquet edge modes 
is equal to the number of edge modes of H with quasi- 
energies between —2it/T and 0. The latter number is 
equal to the total Chern number C of the e = —2n/T 
bands of H. In this way, we deduce that n c d gc = C. 
Combining these two results, the identity |C7| follows 
immediately. 
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